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Abstract 

A simple planar model for an orientational ordering of three- 
fold molecules on a triangular lattice modelling a close-packed 
(111) plane of fullerite is considered. The system has 3- 
sublattice ordered ground state which includes 3 different 
molecular orientations. There exist 6 kinds of orientational 
domains, which are related with a permutation or a mirror 
symmetry. Interdomain walls are found to be rather narrow. 

The model molecules have two-well orientational potential 
profiles, which are slightly effected by a presence of a straight 
domain wall. The reason is a stronger correlation between 
neighbour molecules m triangular lattice versus previously 
considered square lattice 

A considerable reduction (up to one order) of orientational 
interwell potential barrier is found in the core regions of es- 
sentially two-dimentional potential defects, such as a three- 
domain boundary or a kink in the domain wall. For ulti- 
mately uncorrelated nearest neighbours the height of the in- 
terwell barrier can be reduced even by a factor of ICF . 

PACS: 61.48.+C, 7 8. 30. N a 

1 Introduction 

An elegant hollow cage structure of the Ceo fuUerene 
molecule has drawn a close attention of scientists because 
of its unique Ih icosahedron symmetry. A nearly spher- 
ical form of the molecule leads to very unusual physical 
properties of sohd Ceo, fullerite. [3 13 El |3| While at the 
room temperature the molecules can be considered to 
be the exact spheres, the low temperature properties of 
fullerite are determined by the deviation of the molecule 
geometry from the spherical one. At these temperatures 
an orientational molecule ordering takes place, which is 
a basic issue for understanding the results of recent He- 
temperature experiments on heat transport, E] hnear 
thermal expansion, [3 |H| and the specific heat [H] of the 
Ceo fullerite. 

The mass centers of the Ceo molecules in fullerite are 
arranged into an fee structure characteristic for close 
packed spheres with isotropic interactions between them. 
At the room temperature the molecules are found to 
be freely rotating. The resulting crystal space group is 



FmZm. 

Upon the lowering of temperature, the fullerite is sub- 
jected to two transitions. At T « 260 K, it undergoes 
the first order phase transition, after which the fee crys- 
tal lattice is divided into four simple cubic sublattices. 
The molecules are now allowed to rotate about one of 
the 10 molecular threefold axes. Other two of the three 
rotational degrees of freedom are frozen. Within each of 
the four sublattices, the allowed molecule rotation axis 
is fixed along one of the four ([111], [111], [111], or [111]) 
threefold cubic axes, so that the crystal space group is 
Pa3. 

It is worth noting that the truncated icosahedron form 
of a Ceo molecule allows for a more symmetric regular 
crystal structure with only one sublatticc and with the 
four of the above mentioned molecule threefold axes ori- 
ented along the threefold crystal axes (usually regarded 
as a standard molecule orientation). But such a struc- 
ture is energetically unfavourable for the anisotropic 
intermolecular interaction. Instead, the observed low- 
energy structure is obtained by a simultaneous counter- 
clockwise 22° rotation of the Ceo molecules from the ini- 
tial standard orientation about their fixed Pai threefold 
axes. 

As a result of such a rotation, each Ceo molecule is 
oriented with one of the negative charged double C=C 
bonds to every one of the six neighbour molecules be- 
longing to the same close-packed (111) plane perpendic- 
ular to the molecular rotation axis. To the other six 
neighbours (belonging to two adjacent (111) planes) the 
molecule is oriented with the positive charged pentagons 
(P). Therefore following a commonly used notation we 
denote this orientation as "P orientation" . For an ideal 
structure with all the molecules having a P orientation, 
every pair of nearest neighbours is characterized with 
a pentagon from one molecule opposing a double bond 
from another molecule. 

On the other hand, a potential profile of the fullerene 
molecule rotating about its fixed threefold axis has an 
additional metastable minimum^ corresponding to 82° 

^We do not consider to be distinct the energy degenerate min- 
ima obtained by 120° rotation about the threefold molecular axis. 



rotation from the standard orientation (and to 60° ro- 
tation from the P orientation). In this minimum, the 
molecule opposes the neighbour molecules from the same 
(111) plane with the double bonds, and the molecules 
from adjacent planes are opposed with hexagons (H ori- 
entation^). The energy difference between the P and H 
minima is about 11 meV («130 K) and the height of the 
potential barrier is 235-280 meV (w2700-3200 K).[Tn|[n] 

At the high enough temperatures, molecules are able 
to jump between the two energy minima due to the pro- 
cesses of a thermal activation. An average P/H ratio 
is given by the Boltzmann distribution law. Just below 
the high temperature phase transition (Tss260 K) a frac- 
tion of the H oriented molecules is close to 0.5, and for 
Tr:!90 K it tends to 0.15. |I3 

For the temperatures below 90 K the situation changes 
drastically. A waiting time for a molecule to obtain a 
sufficient energy for a jump between the P and H orien- 
tations reaches the order of several days (10** — 10^ s) or 
even more. Therefore at some critical temperature (its 
exact value near 90 K depends slightly on the cooling 
conditions) the molecules become frozen in their cur- 
rent orientational minima, and a transition to an orien- 
tational glassy phase takes place. Below this transition a 
fraction of H oriented molecules remains practically un- 
changed and equal to its equilibrium value (about 15%) 
characteristic for the temperature of the glass transition. 
In other words, in the average every 7th molecule has the 
H orientation, and with the very high probability every 
Cgo molecule has at least one misoriented neighbour. 

While the orientational glass structure is believed 
to persist down to the lowest temperatures, some of 
the experimental data obtained at helium temperatures 
can not be explained with the concept of H oriented 
molecules alone. For example, the data on heat con- 
ductivity OIH] show a maximum phonon mean free path 
of about 50 intermolecular spacings, what implies only a 
0.02 fraction of scattering ("wrong") molecules. Besides, 
the negative thermal expansion _X 8 and the linear con- 
tribution to the specific heat of the fuUerite sam- 
ples at helium temperatures are explained with the help 

■^Strictly speaking, the term 'H (or P) configuration' is more 
adequate for describing a mutual orientation of two neighbour- 
ing molecules. Nevertheless, for every chosen pair of neighbouring 
molecules (let us denote them as A and B) with the fixed directions 
of allowed rotation axes, the mutual orientation depends strongly 
only on the rotation angle of one molecule (say, A). The other 
molecule (B) is always (at any angle of its rotation) turned to 
the first one (A) with a double bond. Therefore, the interaction 
energy of the pair weakly depends on the rotation angle of the sec- 
ond molecule. As to the molecule A, it is at any rotations always 
turned to B with a belt of pentagons and hexagons interconnected 
by single bonds. Thus, namely the molecule A of the pair A,B 
is responsible for the mutual orientation. Aside from this, upon 
the rotation of the molecule A from orientation P to orientation 
H this molecule becomes turned with hexagons (instead of pen- 
tagons) to five more its nearest neighbours. At the same time, the 
energy of its interaction with the other 6 nearest neighbours re- 
mains practically unchanged, because the energy depends mainly 
on the orientation of that latter molecules. Basing on the reasons 
mentioned above, we follow the common notations and use the let- 
ters 'P' and 'H' to denote an orientation of a single molecule, while 
keeping in mind those 6 pair orientations for which this molecule 
rotation angle is crucial. 



of the tunnelling (i.e. quantum) transitions of the Cgo 
molecules between nearly degenerate orientational min- 
ima. Such a possibility was firstly assumed in Ref. |14| . 
where all the molecules in a crystal were supposed to be 
in tunnelling states. However, the paper [13 accurately 
estimates the tunnelling frequency to be about 5.5 K, 
and the number of tunnelling degrees of freedom to be 
~ 4.8 X 10~''(A^/60), where N is the number of carbon 
atoms in a crystal. Obviously, the number of the H ori- 
ented molecules is much bigger, and the above mentioned 
potential barrier between the H and P orientations is too 
high to provide such a low tunnelling frequency. There- 
fore the defect states other than the simple H oriented 
molecules should be considered. 

One of the possibilities for a realization of the low 
potential barrier for Ceo molecule is indicated in our 
previous paper. [T5| Relatively low barrier sites can ap- 
pear within the orientational domain walls because of the 
superposition of the mutually compensating potential 
curves due to interaction with the neighbour molecules 
belonging to different domains. For the case of orien- 
tational ordering of hexagons on a square planar lattice 
considered in jl5j , the height of the potential barrier in 
the wall is found to be 5 times less than in the regu- 
larly ordered lattice. Such a lowering seems to be insuf- 
ficient to provide the necessary magnitude of tunnelling 
frequency following from the available experimental data 
analysis. 231 

Meanwhile, most of the results obtained for a square 
lattice seem to be caused by the incompatibility of the 
molecule threefold C3 symmetry axis with the lattice 
fourfold C4 symmetry axis. In the case of fullerite, a 
fuUerene molecule holds 4 threefold axes and 3 twofold 
axes intrinsic to the fee lattice. Furthermore, the closest- 
packed (111) plane of the Pa3 lattice has a hexagonal 
structure. Six of the 12 molecule nearest neighbours be- 
long to a hexagonal plane, while only 4 of them belong 
to the same square (001) plane. 

Therefore it is interesting and necessary to investigate 
the main features of orientational ordering for the case of 
the molecule symmetry identical to that of the lattice. 
In the present paper we are concerned with the possi- 
ble orientational domain structures formed by the sim- 
ple flat hexagon-shaped molecules arranged into a more 
relevant to a real situation hexagonal lattice, with both 
the molecule and the lattice symmetry axes being C3. 
The main purpose of the paper is to estimate the ener- 
getic inter-molecular interaction barriers for both regu- 
lar close-packed planar structure and for the vicinity of 
extended orientational defects. 

It is a pleasure and a honour for us to devote this paper 
to the Ukrainian low temperature experimentalist of a 
world-wide reputation Prof. Vadim G. Manzhelii whose 
contribution to the physics of cryocrystals in general and 
to the physics of fuUerites and fuUerides in particular is 
well-known and can not be overestimated. 



2 Model 

Let us consider a system of flat hexagonal molecules^ 
located in the sites of a rigid hexagonal planar lattice, 
modelling a (111) plane of the 3D fee lattice. 

Following the empirical potentials used for modelling 
the intermolecular fuUerene interaction (see, for exam- 
ple, Ref. (16j and references therein), we suppose two 
kinds of negative charges, —(1 ± a), to be located at 
the centers of hexagon sides (see large and small filled 
circles in Fig. ^a)). These negative charges recall sin- 
gle and double covalent bonds between carbon atoms at 
the hexagon edge of the truncated icosahedral fuUerene 
molecule. Introduction of the charge parameter a re- 
duces the Cq hexagon symmetry down to C3 intrinsic 
for real Ceo- A requirement of electro-neutrality of the 
model hexagon molecule stipulates a presence of unit 
positive charges at its vertices (shown with the open cir- 
cles in Fig. n^a)). 

For an initial orientation (an analogue of the stan- 
dard orientation in fuUerite) the molecules are chosen 
to be aligned with the positive charges along the lat- 
tice directions. The topmost (positive Y direction) neg- 
ative charge has to be a larger one (see Fig.^a)). The 
molecule rotation angle (j) is measured starting from the 
positive X direction. 

Interaction between the two nearest molecules is given 
with the superposition of the Coulomb interactions be- 
tween the vertices and bonds of these molecules. The 
exact form of the resulting interaction function can be 
found in Ref. (Eqs.(l-4)). Interaction is multipolar, 
so it depends not on the difference of the molecules ro- 
tation angles only (as for the case of the spin systems 
with Heisenberg exchange coupling), but essentially on 
both the angles. So, the energy of interaction of the two 
neighbour molecules characterized with rotation angles 
</>! and (j)2 has to be written as W (j)2) ■ The clock- 
wise and counterclockwise rotations have different effect 
on the interaction: 

l¥(0i; 02) ^ T^(-0i; 02) ^ VF(0i; -02). (1) 

On the other hand, a clockwise rotation of the first 
molecule is somewhat equivalent to a counterclockwise 
rotation of the second molecule. Hence, a combination 
of the lattice mirror symmetry with the molecule mirror 
symmetry leads to the following symmetry relation for 
the interaction function: 

1^(01; 02) = T^(-02;-0i). (2) 

Rotating the molecules shown in Fig. ^a) by an an- 
gle 2tt/3 (or Att/3) about a threefold axis located in the 
center of the triangle 123, one can find that the pair in- 
teraction of the molecules 2,3 (or 3,1) is given with the 
same function W{(f)2; 03) (or W^(03; 0i), respectively). A 

^In some sense, they can be regarded as imitating the Ceo 
molecules viewed along C3 axis. Strictly speaking, such imita- 
tion is competitive only for the fuUerene molecules with the fixed 
C3 axis perpendicular to the considered (111) plane. The Ceo 
molecules belonging to other three Pa3 sublattices have their fixed 
threefold axes tilted to this plane. 
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Figure 1: (a) A local geometry of the model molecules 
on the triangular lattice. Note that molecular rotation 
angles (shown with the help of dashed lines) can be mea- 
sured from any of the three lattice directions, (b) A ge- 
ometry of the simulation cell. Arrows give the 1 — > 2 
order of the input parameters for the pair interaction 
function M^(0i;02). 



relative displacement (which was vertical or horizontal in 
the case of the square lattice considered in Ref. 15 ) of 
the two molecules does not have to be taken into account, 
but an order of the angle parameters is essential. 

For a simulation of the possible domain structures, 
we consider a finite parallelepiped-shaped system, which 
geometry is shown in Fig. ^h). It consists of 20 x 20 
hexagon molecules labeled with the two indexes I and 
m. Arrows show the order of the interaction function 
arguments for each pair of hexagons. The system Hamil- 
tonian then reads: 

N-1 

^ = i^i^'lm.-Al+l.m) 

;,m=o 

,m+l 7 0/m)] 

N-1 

+ YW{<j)iN;(l)l+l,N) 
1=0 



N-1 

+ W{(j)N,m+i;(pNm), (3) 

m=0 

where N — 19, and the last two terms are introduced 
to take into account the edge molecules. For numerical 
simulations, the charge parameter a is chosen to be 1/3. 
A hexagon side makes 0.3 of the lattice spacing. 

3 Possible ordering types 

For a general case of an orientational ordering of the 
identical molecules on a planar hexagonal lattice, the 
structures with 1, 3, 4, and 7 sublattices are possible. 
One sublattice structure would correspond to a uniform 
rotation of all the molecules on a lattice. Three sublat- 
tice structure is characteristic of antiferromagnetic sys- 
tems (Loktev structure 17 ). A close-packed (111) plane 
of the Pa3 structure should contain the molecules be- 
longing to four different sublattices. The results of STM 
imaging of the (111) fuUerite surface JH| confirm this 
fact.'' A more complicated case of seven sublattices could 
be expected for less symmetric molecules. 

As to considered C3 symmetric hexagons, a rather 
aesthetic expectation of the threefold site symmetry^ 
implies either 1 or 3 sublattice case, or a 4 sublattice 
structure involving three identical rotations. Numerical 
calculations give for a ground state a three sublattice 
structure shown in Fig. Efa). It is interesting to notice 
that in the 3 sublattice structure the energy minimum 
corresponds to the molecule positions which do nat pro- 
vide the minimum of the pair potential. The obtained 
molecule rotation angles are {a — 1/3): 

</>! = 72.37209°; 02 = 25.24477°; 03 = 10.87533°. (4) 

Another possible (energy degenerate) ground state can 
be found with the help of the symmetry relation ||2Jl. The 
corresponding angles are given with 

(t>i = -'/'i; 02 = ~<p2; 03 = (5) 

This ground state is related by the mirror symmetry to 
the state defined with Eq. Q). 

The high symmetry of the hexagon molecules makes 
it difficult to perceive the ordering pattern presented 
in Fig. 13 a). A more complicated task of finding an 
orientational defect in this pattern becomes unfeasible. 
Therefore for the purpose of visualization, we implement 
a vector representation of hexagon molecules shown in 
Fig. Elb). The vector rotation angle is three times a 
hexagon rotation angle: = 30;'^. A vector can be 
rotated from 0° to 360°. As a result, a difference between 
sublattices appears to be more clear. 

*At the beginning of the fuUerene era, there were some pubU- 
cations 1191 12UI reporting an 8-sublattice fee structure for the low 
temperature fuUerite. This structure could be obtained by divi- 
sion of each of the four se Pa3 sublattices into two fee sublattices 
with different (P and H) molecular orientations. However, the 
8-sublattice structure has not been confirmed by further investi- 
gations. Therefore, we do not consider it here. 

^An absence of the site symmetry would induce a distortion of 
the lattice. 




Figure 2: A ground state orientational ordering of the 
hexagon molecules (a) , and the same ordering patterned 
with vectors (b). 

Figure |3I shows a change of the molecule interaction 
energy with a change of its orientation for the three 
molecules belonging to three different sublattices. It is 
clearly seen that the molecules are not identical. All of 
them have double-well energy profiles, but the height of 
the interwell energy barrier varies by a factor of 2. The 
potential minima of the 2nd molecule are almost energy 
degenerate, the energy difference being only 1/20 of the 
barrier height (situation, similar to the case of fuUerite). 

For a comparison, the same potential profiles are 
shown for the molecules from the edge of the simulated 
lattice (see Fig. 0J . Such molecules keep only 4 of the 
6 nearest neighbours (molecules from two different sub- 
lattices are missing). As a result, the overall potential 
profile is lowered by a factor of 6/4. The Cgo molecule at 
the fuUerite (111) edge surface is missing 3 neighbours 
from 3 different sublattices. Therefore one could expect 
lowering the orientational barriers by a factor of 12/9.® 

^Nevertheless, it should be emphasized that a change (or rela- 
tively weak lowering) of inter-molecular rotational barriers appears 
to be too small for all the cases of the regular structure to allow 
for an orientational tunnelling which is necessary for a number of 
physical phenomena. One has to remember that the mass of the 
Ceo fuUerene molecule is 720 a.u. It makes very strong constraint 
for the height and the width of energetic barriers which are able to 
give an observable probability (or frequencies 1131 ') of orientational 
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Figure 3: Orientational potential profiles for regular 
molecules belonging to three different sublattices. 
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Figure 4: Orientational potential profiles for the edge 
molecules belonging to different sublattices. 



But the real situation is even more complicated. 
A three-dimensional character of fuUerite lattice leads 
to subdivision of the neighbours of an arbitrary bulk 
fullerene molecule into only two categories, denoted here 
as double-bond (to which the molecule is oriented with 
the double bond), and pentagon (to which the molecule 
is oriented with the pentagon or hexagon) neighbours. 
The six double-bond neighbours belong to the (111) 
plane normal to the molecule fixed C3 axis. The rest 
six pentagon neighbours, which give a major contribu- 
tion to the molecule orientational profile, are located in 
other (111) planes. 

Therefore an edge molecule with a fixed C3 axis nor- 
mal to the edge surface misses three pentagon neigh- 
bours, while the molecules with three other directions of 
allowed rotation axis are missing two double-bond and 
one pentagon neighbour each. As a result, the potential 
relief of a molecule with a normal rotation axis is more 
shallow than the relief of other molecules. In this way, 
the molecules from the four different sublattices which 
are identical by their rotational properties in the bulk 
fuUerite become non-identical at the edge surface crystal 



tunnelling transitions. 



defect due to a loss of the symmetry. This non-identity 
evidently reveals itself in the presence of two additional 
lower temperature order-disorder phase transitions re- 
ported in |21| . 

4 Linear orientational defects 

A general kind three sublattice two-dimensional triangu- 
lar lattice allows for orientational ordering of three dif- 
ferent types. Molecular orientations for these ordering 
types are related to each other with cyclic permutations 
of the rotation angles 0i (« = 1,2,3, cf. Eq. Q) for the 
molecules located at the vertices of a lattice triangle (eg., 
a triangle 123 shown with a solid line in Fig.[2Jl. In the 
case of the considered hexagon molecules an existence of 
the mirror orientational twin defined with Eq. [3 leads to 
appearance of three additional ordering types, which are 
related to the basic permutation ones with the mirror 
symmetry. 

As a result, the considered lattice allows for simulta- 
neous existence of orientational domains with 6 different 
ordering types. A boundary between two domains con- 
tains a linear orientational defect (domain wall). Such 
defect can involve a permutation (clockwise or counter- 
clockwise) or a mirror transformation (with a center at 
3 different lattice sites) of molecular orientations.^ 

A domain wall of the permutation type is presented in 
Fig.E^a). The rotation angles of the molecules located at 
the vertices of a lattice triangle (shown with solid lines) 
have the values (j)2i and (ps in the left domain. In 
the right domain they are equal to </)2, (^3, and 4>i, corre- 
spondingly. The domain wall (grey) is relatively narrow. 
Its width (measured along the horizontal close packed I 
direction) is about one period of 3-sublattice structure. 
As seen along the close packed m direction, this defect 
can be regarded as obtained by removal of one element 
from an ideal sequence . . . 1231231 ... of molecular ori- 
entations. The resulting sequence is . . . 123|231 .... 

Orientational dependence of the potential energy of 
the four central molecules from the domain wall is given 
in Fig. El^b). The molecules are marked in Fig.l^a) and 
labeled with their m index, while I is taken to be 10. 
Molecules 8 and 10 have the orientations of the type 2, 
and the rotation angles of molecules 9 and 11 are close 
to (/)3. The potential profiles are quite close in the form 
to the profiles of the regular molecules (shown in Fig.|3I), 
but one of the two potential barriers is lowered for each 
molecule. 

Orientational domain walls of a mirror nature are 
wider than the permutation ones. Fig. EJa) gives an 
example of the mirror domain wall. For this wall, a se- 
quence of molecular orientations in the m direction is 
. . . 12373 2 1 . . ., where a question mark stands for a 
molecule in the mirror plane. This molecule does not 
fit any regular orientation. Instead, it reflects the mirror 
symmetry of the wall. Figure|Blb) clearly shows the ori- 
entational potential minimum of the molecule 10 to be 

^For the case of fuUerite, there are 4+4 different ordering types 
and 3-1-4 different inter-domain boundaries (not related with the 
symmetry operations). 
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Figure 5: Permutation domain wall (a) perpendicular to Figure 6: Mirror domain wall (a) perpendicular to close- 
close-packed row, and orientational potential profiles (b) packed row, and the potential profiles (b) for the five 
for the four marked molecules identified with the lattice marked molecules with /=10 and with the indicated m 



index m (1=10). 



value. 



located at the rotation angle </> = 60°. Such an orienta- 
tion corresponds to aligning one of the mirror planes of 
the hexagon molecule to the domain wall mirror plane. 

The mirror symmetry of orientational defect is also 
manifested through the symmetry of potential curves of 
other molecules. The potential profiles of the molecules 
9 and 11 (orientations 3 and 3 ), 8 and 12 (orientations 
2 and 2 ) are related through AEg{(l)) — A£'ii(— 0) and 
AEs{(l)) = A£'i2(-0), respectively. 

In the vector pattern of Fig. EJ^a), this symmetry is 
given with the clockwise-counterclockwise vector rota- 
tions on the two different wall sides. Since the rotation 
angles are measured form the X direction, the vectors 
representing molecule 9 and 11 (8 and 12) rotations are 
related with the mirror plane parallel to the X direction. 

The effect of the domain wall on the potential relief of 
the molecules 8 and 12 (orientations 2 and 2 ) is found to 
consist in a slight lowering of one of the two barriers. For 
the molecules 9 and 11 (orientations 3 and 3 ) close to 
the center of the domain wall, both the potential barriers 
are lowered considerably. But the position and the rel- 
ative height of the secondary minimum are unchanged, 
resulting in a shallow character of this minimum seen in 
Fig.EIb). 

The domain walls given in Figs. El have their di- 



rections parallel to one of the sublattice period vectors, 
and perpendicular to one of the close-packed molecular 
row directions. At the same time, there is a possibility 
for a domain wall to lie along the close-packed molec- 
ular rows. An example of a permutation domain wall 
of this kind is presented in Fig. Ufa). The relationship 
between the left and right domains here is the same as 
in Fig. El but the location of the domain wall line is dif- 
ferent. As a result, the molecular orientation sequence 
along the m molecular row can now read . . . 231|312 . . . 
{I = 8), .. .312|123. .. {I = 9), or .. .123|231.. . {I = 10). 
Therefore, the central part of the domain wall contains 
molecules with 6 different potential profiles (orientations 
1, 2, and 3 from the left domain, and orientations 3, 1, 
2 from the right domain). The three potential profiles 
with the lowest energy barriers are shown in Fig. \7\h) . 
Noteworthy that here we gain a low barrier profile with 
almost energy degenerate minima (see dotted curve). 

A mirror nature domain wall parallel to close-packed 
molecular row has a more complicated structure shown 
in Fig.jSfa). It is again wider than the permutation wall, 
so that the molecules from three close-packed rows have 
substantially corrupted orientational potential relief. As 
a result, the number of intra- wall molecules with differ- 
ent orientational profiles increases up to 9, opposed to 
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Figure 7: Permutation domain wall parallel to close- 
packed row (a), and orientational potential profiles (b) 
for the marked molecules with the given lattice indexes 
/ and m. 

6 different profiles for a permutation wall. Furthermore, 
the direction of the domain wall does not coincide with 
the lattice mirror plane, so there are no mirror symme- 
try in the pattern of Fig.lHl^a), and, respectively, no sym- 
merty relations for the potential curves (conf. the mirror 
symmetry of the potential profiles shown in Fig. Elb) for 
the domain wall perpendicular to close-packed molecu- 
lar row). In Fig.lHIb) we give the orientational potential 
profiles for the three molecules with the lowest interwell 
energy barriers. It is seen that there exists a molecule (o) 
which interwell barrier is about 1.4 times lower than the 
lowest of the regular molecules interwell barriers. The 
molecule is situated in the center of the domain wall 
and marked with a circle. The corresponding potential 
profile is plotted with a sohd fine in Fig.lHl^b). The ob- 
tained reduction of the orientational interwell barrier is 
caused by a less correlation between the nearest neigh- 
bour molecules (every molecule within the considered 
wall has neighbours of six different orientations). 

5 Two-dimensional defects 

The results on the modelling of the straight domain walls 
in the considered system show that the molecules with 
the most shallow potential profile tend to appear at sites 
with the reduced correlation between the orientations of 
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Figure 8: Mirror domain wall parallel to close-packed 
row (a), and the potential profiles (b) for the marked 
molecules (lattice indexes I and m are indicated). 

the neighbour molecules. For the straight walls such 
condition is met at the boundary of two domains with 
different sets of equilibrium molecule orientations (mir- 
ror domain walls). 

Then it is straightforward to continue the search for 
the shallow potential molecules in the core regions of 
essentially two-dimensional orientational defects. One 
of such promising two-dimensional defects is a meeting 
point of three different domains. Molecules at the cen- 
ter of this defect should have three pairs of neighbours 
belonging to three different domains, so one could ex- 
pect for an additional decrease of interwell barriers with 
respect to the two-domain boundary case. 

The results of numerical calculations indeed show the 
further reduction of interwell potential barriers at the 
boundary of three orientational domains. The most ef- 
fective reduction is found to take place in the presence 
of mirror boundaries. 

Figure ISJ^a) shows an example of orientational defect 
formed at the intersection of three domain walls per- 
pendicular to molecular rows. The left (narrow) domain 
wall is of a permutation type, the other two (the bot- 
tom one and the right one) have a mirror nature and are 
much wider. The right domain wall incorporates a kink 
in order to minimize a surface spanned by the defect. 
Molecules with the lowest interwell barriers are marked. 

As it is mentioned above, no significant potential 
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Figure 9: Structure of orientational defect (a) formed 
at the boundary of three orientational domains, and the 
potential profiles (b) for some chosen molecules in the 
defect core region. Pairs of molecules marked with the 
same sign (A or O) have the symmetry related potential 
profiles, so only one of the profiles is given for each pair. 
Molecules are labeled with I and m indexes. 



barrier reduction has been observed for straight do- 
main walls perpendicular to close-packed molecular rows. 
Therefore the marked molecules can be seen only at the 
crossing of the three walls. The corresponding orienta- 
tional potential profiles are plotted in Fig. I^b). 

It is surprising that the potential profile with the least 
energy barrier belongs not to the molecule □ situated 
in the very center of the defect (potential curve plotted 
with a dashed line) with totally different orientations 
of all the 6 nearest neighbours, but to the molecule lo- 
cated at the beginning of the bottom domain wall (o, 
solid line). For the last molecule the orientations of the 
nearest neighbours differ only slightly from that in the 
straight wall, but the interwell potential barrier is 2.3 
times lower than the lowest regular molecule barrier. 

The other four molecules which are marked in Fig.lHIa) 
are located within the center of the kink in the right 
domain wall. At a closer look, one can find a kind of 
a symmetry center at the middle of the line between 
the molecules marked with A. An exact symmetry is 
following: if the centers of two molecules are related 



Figure 10: The structure of the kink in the mirror do- 
main wall (a), and orientational potential profiles (b) 
for marked molecules. Only one potential curve is given 
for every pair of symmetry related molecules which are 
marked with identical signs. Indexes 1 and m are indi- 
cated. 



with inversion symmetry, these molecules have the ro- 
tation angles which have equal absolute values, but dif- 
ferent signs. Therefore the two molecules marked with 
A (as well as the two molecules marked with O) have the 
same orientational dependence of intermolecular interac- 
tion potential, the only difference being in the clockwise 
or counterclockwise direction of molecule rotation. This 
can be compared to the symmetry of potential curves in 
Fig.EIb), but there is no mirror plane in the present case. 
To avoid having a very complicated picture, only one of 
the two symmetry related curves is shown in Fig.l^Jb) for 
each pair of molecules. Both the dotted and the dash- 
dotted curves have an interwell energy barrier which is 
less than the lowest energy barrier characteristic for reg- 
ular molecules. This means that at the center of the 
kink in a domain wall (also a two-dimensional defect) 
molecules have ill-correlated nearest neighbours. There- 
fore the case of a kinked domain wall has to be investi- 
gated more thoroughly. 

Figure ITUT a) shows a structure of the kink which con- 
tains the molecule with the lowest height of the orienta- 
tional interwell barrier obtained in our simulations. This 



molecule (in fact, two molecules, since the kink has a cen- 
ter of symmetry of the described above kind) is located 
at the very center of the kink, and the corresponding 
potential curve is shown in Fig. llOf b') with a solid line. 
The height of the interwell potential barrier is already 5 
times less than for the case of regular molecules. 



6 Totally uncorrelated neigh- 
bourhood configuration 

The three-dimensional defect structure of the real ful- 
lerite can be even more complicated. As a result, some 
molecule can have the neighbours which orientations are 
fixed by different elements of the defect network. In 
the frame of our simple two-dimensional model such 
neighbourhood would be totally uncorrelated, and the 
height of the interwell barriers could be further lowered. 
Therefore it is interesting to know a minimum possi- 
ble height of the molecule interwell potential barrier for 
an arbitrary orientational configuration of its neighbour 
molecules. 

For this purpose, let us consider a system of 7 hexagon 
molecules located at the sites of hexagonal lattice, so 
that one central molecule has 6 nearest neighbours. Ro- 
tation angles of the outer molecules are fixed to be equal 
to 6 random numbers between 0° and 120°, and then the 
orientational potential profile of the central molecule is 
calculated. Configurations with the most shallow poten- 
tial profiles obtained in the course of about 10® different 
realizations of random neighbourhood configuration are 
shown in Figs. [Til 113 andlHl 

Figure 1111 gives an example of a molecular config- 
uration with interwell potential barriers of a central 
molecule reduced by two orders of magnitude with re- 
spect to the case of totally orientationally ordered lat- 
tice. This configuration is nearly symmetric (the outer 
molecules have rotation angles about ±30°). The cen- 
tral molecule has a four-well orientational potential pro- 
file with the main minimum located close to 30°. One 
could expect that a completely symmetric configuration 
might have even more shallow potential profile of the 
central molecule, because of the increase of the inter- 
action energy at the minima of the potential. Contrary 
to the expectations, the exactly symmetric configuration 
(not shown) has an order of magnitude higher interwell 
barriers than the one shown in Fig. 1111 Thus, interwell 
barriers prove to be extremely sensitive to even very smal 
rotations of the molecules. 

The case of the molecular configuration with a two- 
well orientational profile of the central molecule is pre- 
sented in Fig. lEl If one does not take into account the 
difference between the values of negative charges, this 
configuration seems to be close to having a mirror sym- 
metry. Probably, namely this difference leads to an in- 
crease of the interaction energy at the potential minima. 

The molecular configuration with the lowest obtained 
interwell potential barrier of the central molecule (shown 
in Fig. I13|l has no symmetry at all. The orientational 
profile has three minima of different depth, while the 
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Figure 11: A molecular configuration with nearly sym- 
metric orientations of the outer molecules (a) and the 
corresponding shallow potential profile of the central 
molecule (b). An inset in the bottom panel shows a 
magnified potential curve. 

lowest interwell barrier is about 200 times lower than 
the corresponding lowest barrier in the regularly ordered 
lattice. 

Also it should be noted that the molecules of the reg- 
ularly ordered lattice (namely, the molecules with the 
(pi orientation, see Fig. |2Jl have the neighbourhood con- 
figuration with the highest possible interwell potential 
barrier. While minimizing the overall interaction energy, 
this configuration minimizes also an interaction energy 
at the minimum of the one molecule potential, and deep- 
ens this minimum. 



7 Discussion 

The considered simple planar model recovers some of the 
features of the fuUerite lattice. First of all, it predicts a 
multi-sublattice structure for a system which would be 
arranged into a more symmetric 1-sublattice manner in 
the absence of anisotropic intermolecular interactions. 

Then, the model involves lowering of orientation po- 
tential relief of the molecule at the crystal surface. This 
can be compared favourably to the absence of H oriented 
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metric configuration with a two-well potential profile. 



figuration with the lowest found interwell barriers. 



molecules at the STM image of the fullerite surface. [TH| 
Futhermore, at a closer look this image shows a slight 
difference in orientations of the fullerite molecules be- 
longing to the three sublattices which have the molecular 
C3 rotation axes not perpendicular to the surface. This 
difference is due to a competition of the two-dimensional 
character of a surface (with probably another subdivision 
into sublattices) with the bulk equilibrium orientaitons 
of the molecules below the surface molecular layer. 

Rather narrow character of the domain walls in the 
considered high-symmetry system is rather natural for a 
system with only one kind of interaction involved. * It 
agrees well with a very sharp character of a domain wall 
observed in a two-dimensional monolayer of Cgo fuUerene 
molecules. [22 This wall contains also a very sharp kink, 
which is a kind of essentially two-dimensional defect that 
can incorporate molecules with low orientational inter- 
well barrierrs. 

The sharp character of the observed kink implies a 
possibility of existence of strongly localised orientational 
defects also in the bulk of the three-dimensional fullerite. 



*For the case of ferromagnets the domain wall width is of the 
order of ay/ J/A, where a is a lattice spacing, J is an exchange, 
and A is an anisotropy. Since A is a relativistic correction, the 
ratio J/A can be increased up to 10®. But for the present case of 
one interaction this ratio is about 1. 



Some of this strongly localised defects with necessarily 
uncorrelated orientations of the neighbour molecules ori- 
entations should involve molecules with an orientational 
potential which is sufficiently shallow to give a reason- 
able frequency of tunnelling transitions. As to the rather 
high Cgo molecule mass, the recent molecular dynamics 
simulations on the dislocation kink tunnelling '23' in Ag 
show an efficiency of tunnelling of complex heavy object 
under certain conditions. 

The idea to explain negative thermal expansion of 
solids with the double-well tunnelling statistics was 
suggested by Freiman in 1983 for the case of solid 
methane. |24| In the range of temperatures where the 
conventional phonon mechanism does not work the ther- 
mal expansion is established as a result of competition 
of two factors. The first factor is a lattice contraction 
due to the process of populating the tunnel states with 
an increase of temperature. Shrinking the distances be- 
tween molecules increases the height of the orientational 
interwell barriers, what leads to a decrease of the tun- 
nelling energy splitting, and as a result to a decrease of 
the system free energy. The contraction of the lattice 
is stabilised by an increase of an elastic part of the free 
energy for every fixed value of crystal temperature. 

Since the population of tunnel states has very strong 
exponential temperature dependence, the thermal ex- 



pansion resulting from the competition of the two factors 
is always negative. At T w K it is practically absent 
(no molecules on excited tunnel levels) . With an increase 
of temperature the population of the excited state grows, 
therefore the lattice is contracted. But at T > A, where 
A is the tunnel state energy splitting, both the ground 
and the excited states become almost equally populated, 
so that the effect becomes much less pronounced. This 
means that there should exist a maximum in the magni- 
tude of the negative thermal expansion coefficient. 

With a simple differentiation of expression (6) of 
Ref. |24| one can obtain that this maximum takes place 
at the temperature T^ax satisfying the equation 

2T„,ax/A - tanh (A/2T^ax + 1/2 In (/1//2)), 

where /i and /2 are the degeneracies of the ground and 
the excited states, respectively. It is easy to see that 
Tmax < A/2 holds for any ratio /1//2. 

Therefore the tunnelling energy splitting in fuUerite 
can be estimated from the Tmax position in Refs. [71 |S| 
to be more than 8 K. On the other hand, the positive 
thermal expansion of pure fuUerite at T < 2 K implies 
a presence of processes other than two-well tunnelling 
(probably, the conventional phonon mechanism is still 
valid) at this low temperature. 

The possibility to detect experimentally the negative 
contribution to thermal expansion due to the tunnelling 
objects depends strongly on the relative magnitude of 
the positive (conventional) and negative (tunelling in 
this case) contributions. In the case of fuUerite, the 
negative contribution is more pronounced, but one still 
encounters a difficulty in determining the tunnelling ob- 
ject. The first hypothesis about a tunnelling of regular 
Ceo molecules between P and H orientations IT had 
a drawback of high interwell potential barrier. Further 
introduction of the idea of competition of isotropic and 
anisotropic parts of intermolecular interaction potential 
(though in orientational glass) 25 has lead to a current 
understanding (given in our previous paper \T5\ and the 
present one) that the tunnelling objects are to be found 
at the strongly localised orientational defects of the ful- 
lerite structure. The more detailed description of such 
defects could be obtained with the help of more realistic 
three-dimensional modelling of the Ceo crystal structure, 
what should be a subject for future studies. 
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